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I. Introduction 



Ecology is the study of inter-relationships between organisms and environment. It is natural that two 
or more species living in a common habitat interact in different ways. The application of mathematics to 
ecology dates back to the book “An Essay on the Principle of Population” by Malthus. In recent decades, many 
researchers [1,2] have focused on the ecological models with three and more species to understand complex 
dynamical behaviors of ecological systems in the real world. They have demonstrated very complex dynamic 
phenomena of those models, including cycles, periodic doubling and chaos [6,7]. The discrete time models 
which are usually described by difference equations produce much richer patterns [3, 4, 5]. Discrete time 
models are ideally suited to describe the population dynamics of species, which are characterized by discrete 
generations. 



II. Mathematical Model 

In this paper, we consider the one prey- two predator systems describing the interactions among three 
species by the following system of difference equations: 



J („ + i) = r JW (i-^))- ‘ K( " )j,( " ) 

1 + x(n) 



bx(n)z(n) 
1 + x(n) 



y(n + T) = y(n)(l-c) + 



ax(n)y(n) 
1 + x(n) 



( 1 ) 



z(n + 1) = z(n)(l- d) + 



bx(n)z(n) 
1 + x(n) 



where x(n) Prey, y(n) and z{n) are two predator respectively, and all parameters are positive constants. 
The parameter a is the intrinsic growth rates of the prey population, b denotes the death rate of the mid- 
predator, c denotes the death rate of the top predator, d is the rate of conversion of a consumed prey to a 
mid-predator, e is the rate of conversion of a consumed mid-predator to a top predator. 



III. Existence Of Equilibrium 

The equilibrium points of (1) are the solution of the equations 
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The equilibrium points are E 0 =(0,0,0) , E x - 



1 






i--,o,o 

V r 



. E 2 = 



d d(\-2r)+ (1+ r)b 



-d + b 



(b-d) z 



IV. Dynamical Behavior Of The Model 

In this section, we study the local behavior of the system (1) around the equilibrium points. The 
stability of the system (1) is carried out by computing the Jacobian matrix corresponding to each equilibrium 
point. The Jacobian matrix J for the system (1) is 



J(x,y,z) 



on ay + bz axy + bxz 
r( 1 - 2x) - — + - 



1 + x 



a+*r 



-ax 
1 + x 



ay axy 



1 + x (1 + x) 2 
bz bxz 



1 -c + 



ax 
1 + x 



-bx 
1 + x 

0 



0 



1 — d + 



bx 



(2) 



1 + x (1 + x) 1 + x 

Theoreml. The equilibrium point E 0 is locally asymptotically stable if0<r<l, 0<c <2 and 0 < d < 2 
otherwise unstable equilibrium point. 

Proof: The equilibrium point (0,0,0) in the system (1) is asymptotically stable. The Jacobian matrix for E () 
is evaluated as 



J(E 0 ) = 



r 0 0 

0 1-c 0 

0 0 1-d 



The Eigen values are \ 2 3 = {r,l — c,l — r/} 

If the equilibrium point E 0 satisfying the existence condition then it is most have one the following possible 
types: 

(1) Sink if 0<r<l , 0<c<2 andO <d <2 . 

(2) Source if r > 1 , c > 2 and d > 2 . 

(3) Non-hyperbolic if r — 0 , c — 2 and d = 2 . 

(4) Saddle otherwise. 



Theorem 2. The equilibrium point E x is locally asymptotically stable ifl<r<3, 



a{r- 1) 
2r — 1 



<c< 2 + 



a(r-l) 
2r — l 



and 



b(r- 1) 
2r-l 



<d< 2 + 



b(r- 1) 
2r-l 



otherwise unstable equilibrium point. 
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Proof: 

Corresponding to the equilibrium point E x , the Jacobian is 

a(r- 1) 



J(E X ) = 



2-r 



0 l-c + 



2r-l 

a(r-l) 



2r — l 



1 — d + 



b(r-l) 

2r — l 

0 

bjr — 1) 
2r — l 



m „ air — 1) „ , b(r-l) 

The eigen values are A = 2-r,l-cn ,1 -d-\ 

%2 ’ 3 2r — 1 2r — 1 

If the equilibrium point satisfying the existence condition then it is most have one the following possible 
types: 

(1) Sink if l<r<3, 

(2) Source if r > 3 , 



a(r-l) „ a(r-l) 

<c< 2 + 


and*"” 


^ Hr- 1) 

<d < 2h . 


2r — \ 


2r — 1 


2r — 1 


2r — 1 


a(r-Y) „ a(r — 1) 

> c > 2 + 


and Hr ~ l) 


, , „ *(r-l) 

> d ox d > 2 H 


2r — 1 


2r-l 


2r — l 


2r-l 



(3) Non-hyperbolic if d = 1 and c = orc , 2+ «L± and d or rf -2+*^ . 

2r — 1 2r — 1 2r-l 2r-l 

(4) Saddle otherwise. 

Theorem 3. The equilibrium point E 2 is locally asymptotically stable if and only if A , C , and AB - C are 



positive. 

Proof: The Jacobian matrix J at E 2 has of the form 



J(E 2 ) = 



r(b-2d)-rd f J(-2r + l) + (l + r)Z?^ 



b-d 
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-ad 
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b ’ 



( 3 ) 



„ . ad 

a 2X — (J , a 22 — 1 — c H , a 23 

b 



= 0 , a 31 



d(-2r + l) + (l + r)b 

b 



&22 ^ 9 ^33 1 



Its characteristic equation is 
A 3 + AA 2 + BA + C = 0 (4) 

with A=-(a n +a 22 +a 33 ), B = a ll a 33 + a 22 a 33 + a n a 22 - rz 23 a 32 and C=(a 23 a 32 - a 22 a 33 ) a n . By the Routh-Hurwitz 
criterion, E 4 is locally asymptotically stable if and only if A , C , and AB - C are positive. 

V. Numerical Simulations 

Figure 1: The parameter values r- 0.99 , a = 0.95 , Z? = 0.9 , c = 0.97 , J =0.79 , the system (1) has an 
equilibrium point E 0 in the Time Series Plot and Phase Portrait which is globally asymptotically stable. 
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Fig 1. Time Series Plot and Phase Portrait at E 0 



Figure 2: The parameter values r=2.89, a=0.9, b=1.35, c=0.99, d=1.09, the system (1) has a equilibrium point 
Ei in the Time Series Plot and Phase Portrait which is globally asymptotically stable. 




Fig 2. Time Series Plot and Phase Portrait at Ei 

Figure 3: The parameter values r=1.98, a=1.78, b=0.41, c=0.94, d=0.09, the system (1) has a equilibrium point 
E 2 in the Time Series Plot and Phase Portrait which is globally asymptotically stable. 




Fig 3. Time Series Plot and Phase Portrait at E 2 
VI. Bifurcation 

Bifurcation is a change of the dynamical behaviours of the system as its parameters pass through a 
bifurcation (critical) value. Bifurcation usually occurs when the stability of an equilibrium changes. An 
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examination of the bifurcation diagram of the system shows that a given qualitative change in the functioning 
of the ecosystem may be produced by changing different parameters. We restrict our analysis to events that 
occur when the parameter a is changed, and the other parameters are fixed at given values. In this section, we 
focus on exploring the possibility of chaotic behaviour for the prey population respectively. 
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Fig 4. Bifurcation diagram for prey population 



Variation of prey population 
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Fig 5. Variation of prey population 

REFERENCES 



Journal Papers: 

[1] . Abd-Elalim A. Elsadany, Dynamical complexities in a discrete-time food chain, Computational 

Ecology and Software , 2012, 2(2): 124-139. 

[2] . Sophia R. -J. Jang, Jui- Ling Yu, Models of plant quality and larch hudmoth interaction, Nonlinear 

Analysis(2009), doi: 10.1016/.na.2009.02.091. 

[3] . Debasis Mukherjee, Prasenjit Das & Dipak Kesh, Dynamics of a Plant -Herbivore Model with Holling 

Type II Functional Response. 

[4] . A. George Maria Selvam, R. Dhineshbabu, and V.Sathish, Dynamical Behavior in a Discrete Three 

Species Prey -Predator System , International Journal of Innovative Science, Engineering & Technology 
(IJISET), Volume 1, Issue 8, October 2014, ISSN: 2348-7968, pp 335 - 339. 

[5] . A. George Maria Selvam, R. Janagaraj, P. Rathinavel, A Discrete Model of Three Species Prey- 

Predator System, International Journal of Innovative Research in Science , Engineering & Technology, 
ISSN (Online): 2319-8753, ISSN (Print): 2347-6710, Vol. 4, Issue 1, January 2015, DOI: 
10. 15680/IJIRSET. 2015. 0401023, pp: 18576 to 18584. 

[6] . Prasenjit Das, Debasis Mukherjee & Kalyan Das, Chaos in a Prey- Predator Model with infection in 

Predator- A Parameter Domain Analysis , Computational and Mathematical Biology. 

Books: 

[7] . Saber Elaydi, An Introduction to Difference Equations, Third Edition, Springer International Edition , 

First Indian Reprint, 2008. 



I IJMER I ISSN: 2249-6645 I 



www.ijmer.com 



I Vol. 5 I Iss.3l Mar. 2015 I 541 



